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Abstract 

We discuss the ,f{R) gravity model in which the ori- 
gin of dark energy is identified as a modification of grav- 
ity. The Noether symmetry with gauge term is investi- 
gated for the f{R) cosmological modeL By utilization 
of the Noether Gauge Symmetry (NGS) approach, we 
obtain two exact forms f{R) for which such symmetries 
exist. Further it is shown that these forms of f{R) are 
stable. 

Keywords f{R) gravity; cosmology; Noether symme- 
tries. 



1 Introduction 

From th e observational data of supernovae of T ype la 
(SN la) JRiess et allll998l : IPerlmutter et allll999t ). it was 
reported that the present observable Universe is under- 
going an accelerating phase. The mysterious source for 
this late-time acceleration was dubbed 'dark energy'. 
Despite many years of research (see e.g., the reviews 
( Copeland eta3l200d[L"uca fc Tsuiikawall2010l )) its ori- 
gin has not been identified yet. Dark energy possesses 
negative pressure leading the accelerated expansion of 
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the Universe by counteracting the gravitational force. 
One possible approach for the construction of dark en- 
ergy models is the modification of the geometrical part 
of the Einstein equations. This approach known as 
Modified Gravity can successfully explain the rotation 
curves of galaxies, the motion of galaxy clusters, the 
Bullet Cluster, and cosmological observations without 
the us e of dark ma tter or Einstein's cosmological con- 
stant |Moffatll20lil) . 

The f{R) theory of gravity is a candidate of modified 
theories of gravity which is obtained by replacing the 
Ricci scalar R with an arbitrary fu nction f(R) in the 
Einstein- Hi lbert Lagrangian (see (JNoiiri fc Odintsov 
20071 I2OIOI ) for reviews on this topic). Recently dif- 
feren t forms of f(R) have been prop osed rB amba et all 
2008t lElizalde et allboill : ICognola et al 20081 ). These 
theories can produce inflation, mimic behavior of dark 
matter and current cosmic acceleration if suitable 
positive and negative powers of cu rvature are added 
into the Einstein- Hilbert action (INoiiri fc Odintsov 



2003 : lBamba et alll2008l:ICarroU et alll2004l:IStarobinskv 



20071 ). Besides, compatibility with observational data, 
the modified gravity theory has to be compatible 
with general relativity and other viability conditions 



( Noiiri fc Odints ovl2006ll2007tlLuca fc Tsuiika wa'2010l:l 
Setare fc Jamil 201l HTsuiikawal2010l:IJamil et al 20f l|. 
The plan of this paper is as follows: In Section II, 
we discuss briefiy the basics of the f{R) model and con- 
struct the modified field equations. Then in Section III, 
we discuss the Noether gauge symmetry approach for 
the f{R) modified gravity. In Section IV, we check the 
stability conditions on the f{R) functions and finally 
conclude our paper in the last Section. Throughout 
this paper, we choose units c = 1 = IGttG and the 
metric signature (— , -I-, -|~, +). 



2 Field equations in f{R) theory of gravity 

In this section a spatially flat Friednian-Robertson- 
Walkcr (FRW) Universe within the framework of f{R) 
gravity is considered. Our aim is to investigate mod- 
els which exhibit Noether symmetry with gauge term; 
we do not take in account matter contribution in the 
action. We start with a (3+l)-dimensional action 
(iFaraoni fc NadeaullJOOsI) 



S 



d^x. 



-gf{R), 



(1) 



where R is the scalar curvature and f{R) is an arbitrary 
non-linear function of R. The first model proposed had 
the form f{R) = R — M^/-R, in which the correction in 
R^^ becomes important only at low curvatures R — > 
0. It is demonstrated in (jStellelll977l : iBuchbinder et al 
19921 ) that action (1) is never renormalizable. In order 
it to be renormalizable, as it is shown there, L should 
consists of i?^ + R^uR''" + R + A. In addition to the 
desired phenomenological properties of modified gravity 
in cosmology, th ere is some motivation fo r these models 
from M-thcory (JNoiiri fc Odintsovlliooi ). 

Variation of the action (1) with respect to the metric 
yields the field equations 



1 



9,..f{R) - R^ufiR) + V^V,/'(i?) 
-3M-n/'(i?) = 0, 



(2) 



where a prime denotes derivative with respect to R and 
H, V = 0, 1,2,3. The operator V^ represents covariant 
derivative and □ = VV^. We assume that the ge- 
ometry of the spacetime is given by the flat FRW line 
element 



ds^ 



-df + a^{t){dx^ + dy^ + dz^). 



(3) 



With the consideration of this model the field equations 
become 



2H + 3H^ 



H'- 



1 

7 



f'R' + f"{2HR + R) 



+^{f-Rf') 



6/ 



-[f'R-f-GRHfl 



(4) 
(5) 



where H = a/a is the Hubble parameter and overdot 
denotes the derivative with respect to proper time t. To 
investigate NGS an effective Lagrangian for the model 
is needed whose variation with respect to its dynam- 
ical variables yields the correct form of the equations 
of motion. We follow the work of Souza and Faraoni 



( de Souza fc Faraonil 120071 ) and consider the above ac- 
tion that represents a dynamical system in which the 



scale factor a and the scalar curvature R play the role 
of dynamical variables. The actio n (1) can be written 
as (jCapozziello fc Lambiaselbood ) 



S = I dtC{a, a, R,R)^ I dt a^f{R) 



(6) 



where the definition of R is introduced in terms of a 
and its derivatives as a constraint. In order to apply 
the NGS approach, one may easily verify that, in the 
FRW model, the Lagrangian related to the above action 
(6) takes the form 

C{a, a, R, R) = 6{d^af + aRa^ f") + a^{f'R - /). (7) 

Varying the Lagrangian (7) with respect to R yields the 
following relation for the scalar curvature 



i? = 6 



{"'-f)- 



(8) 



The equation (4) can be obtained by varying the La- 
grangian (7) with respect to a. 



3 Noether gauge symmetries in f{R) theory of 
gravity 

Noether symmetries are the symmetries of the La- 
grangians which have found a recent impetus since 
these can reveal new features of the gravitational the- 
orics iCapozzicllo &: Lambiasd 2000t Capozziello et all 
r2009b. 2007: Capozziello fc Lambiasell2007l ). The NS 
are essenti al tools for solvin g the gravitationa l field 
equations (jSanval et all 120051 : IJamil et all 1201 ll ). In 
Scalar- Tensor cosmology, this approach results in an 
extra correction term R^^ and fixes the form of the 



coupling parameter and the field potential (JMotavalia 
20081) . This approach yields an exact for m of fJR) fun c 



tions relevant for cosmological mode l (Vakili 2008allb . 
2010|: ICapozziello fc De Fehc j l2008t ICapozziello et al 



2007 



Roshana fc Shoiail |2008|). In the literature. 



Noether symmetries have been studied in the context of 
f{R) theory of gravity by ignoring the gauge fun ct ion in 
the Noether symmetry condition JVakililboOSalbl. 12010 ; 
ICapozziello fc De Felic^boOSt [Capozziello et alll2007l) . 
We consider the gauge term of the Noeth er symmetries 
( Ibragimovlll999l : iBluman fc Kumeilll989l ) . Here we ap- 
ply the approach of NGS to look at some interesting 
forms of f{R)- 
A vector field 

o o o 

l^ = ^(t,a,R)-+f^{t,a,R)—+p{t,a,R) — , (9) 



whose first prolongation is 

Xl^l = X + r)(i, a,R)—+ ${t, a,R)-^, 

where 



(10) 



. _ dp 80 -dp dc ^2di 

-aR—, 
oa 



(11) 



is called a NGS if the following condition holds 
xW/: + (De)/: = DA(t,a,-R). (12) 

Here A is the gauge function and 

(13) 



T^ 9 ■ d ■ d 



Using the Lagrangian (7) in (12) and after the sepa- 
ration of monomials we obtain the following system of 
determining equations 



ia = 0, in = 0, 



(14) 



77/' + /3a f + 2af'ria - itaf + a'f"p,a = 0, (15) 



12af'iit+6a'f"l3,t^A,a, 



rv,R = 0, 



6a2/"77,t = A,fl, 



(16) 



(17) 



(18) 



2ar7jWr'f3Wf'{v.a+l3,R-it)+2af'Tj,^n = 0, (19) 



which leads to a contradiction. Hence we take /" ^ 0, 
so from ITTI) 



V,R = 0. 



(23) 



There are three cases in the solution of the above de- 
termining equations (14)-(20). These are given below. 

Case-I. If / is arbitrary, the solution of the above 
system (14)-(20) gives rise to the following NGS gener- 
ator 



X- ^ 



(24) 



and the gauge term is a constant which can be taken 
as zero. The energy type first integral is 

/ = Gd^af + QaRa^f" - a^{f'R - J). (25) 

Case-II. If / is the fractional power law, viz. 

/ = kR"\ (26) 

then the solution of the determining equations (14)-(20) 
yields 



■^ —1 

77 = -61a + 63a , 
li = -2R{hi+h^a-'^). 



(27) 



The gauge term is zero. The corresponding generators 
are 

X - ^ 

d 2 d d 

^^ = 'Wt^r^a-^^'mi^ 



(28) 



which constitute a well-known three-dimensional alge- 
bra with commutation relations 



a\iij + aS„t)[fR -f) + a^f'RP = A,t 



(20) 



where "," denotes partial derivative. 

From ([17]), for /" — 0, we get 77. /j ^ 0, then from 
(|19p . /' = 0. Therefore the two functions ^ and A 
become only functions of the variable t and ((20|) reduces 
to 



a2(37? + aG)/-Af 



(21) 



Differentiating this last equation with respect to R we 
get 



'7,H = 0, 



(22) 



[Xo,Xi] = 0, [Xo,X2] = 0, [Xi,X2] = --X2. 

Here too the gauge term is zero. The corresponding 
first integrals are 

h = -9d^atfnR^/^ + -f^ta^R^'^ + iaa^f^R^/'^ 

+3Ra^faR-^/^ - ^aRta^foR^^^^, 
h = 9afoR^/^- + ^afoR~'^^R. 



(29) 



Ca se-Ill. If / is a general power law of the 3 



form dMotavalial boOSt ICaDOzziello fc De Felici I2OO8 : 



CaDOzzicUo & Lambiasc 200(j; Capo zziello et alll2009b . 



2007: Capozzicllo &: LambiasQ;2007i) 



/ = foR\ 



'^y^O,!,-, 



(30) 



then the above linear determining equation system has 
solution 

i = ht + b2, ij = ^^-^abi, p^-2Rbi. (31) 



Here the gauge term is zero and /o is a constant in (30). 
Note that if i^ = 0, 1, |, then / is a constant, linear or 
fractional power law of Casc-II. These are thus excluded 
in Case-Ill. It should also be stated that for Casc-III 
there are two Noethcr symmetries. The Noether sym- 
metry generators are given by 






2R 



d_ 
OR' 



(32) 



The Lie algebra is the two-dimensional Abelian Lie al- 
gebra, [Xo,Xi] = 0. We note that for the fractional 
power law the algebra is three-dimensional and for the 
arbitrary power law it is two-dimensional and thus the 
symmetry breaks by one. 

The corresponding first integrals for this case are 



/o 



/i = 



Qa^avf^R"-^ + QaRa^v{v - l)faR''~^ 

-a^foR^y-l), 

-Ga^atfovR"-^ + ahfo{iy - l)R'' 

+Av{2 - v)aa^fQR''-^ + 2v{v - l)(2i/ - 1) 

xRa^foR''-^ - 6i^{iy - l)ama^ foR"-'^ . (33) 



Hence we have obtained some different f{R) func- 
tions from our analysis. In the next section, we check 
the conditions under which our f{R) functions can pro- 
duce viable cosmology. 



4 Stability Analysis 



f{R)^R-2A for i? » i?o- This is required for the 
local gravity tests and for the presence of the matter 
dominated epoch. 
< -^(r = -2) < 1 at r == ^^ = -2. This is 



required for the late-time de Sitter point. 

To ensure classical and quantum stability in the re- 
gion R, we want our /(i?) theory to satisfy condi- 
tions (1) and (2). The first condition means that grav- 
ity is attractive and the graviton is not a ghost. It 
was recognized long ago that its violation during the 
time evolution of a FRW background results in the im- 
mediate loss of homogeneity and isotropy and forma- 
tion of a strong spac e-like anisotropic curvature sin- 
gularity (|Nariailll973l ). The second condition is nec- 
essary so that in the case of f{R) models of present 
dark energy, the necessity to keep it valid for all val- 
ues of R during the matter- and radiation-dominated 
epo chs in order to avoid the Dolgov-Kawasaki in stabil- 
itv (iNoiiri fc Odintsovl[2003ll2004t liaraonilboOGl) . Con- 
dition (3) imposes that the modified f{R) gravity must 
reduce to Einstein gravity under a suitable limit of the 
curvature scalar. Finally condition (4) tells us of a de 
Sitter point which corresponds to a vacuum solution 
with constant R i.e. D/' = at this point. It triv- 
ially implies that any quadratic f{R) ^ R? will satisfy 
this condition and gives rise an exact de Sitter solution. 
Hence quadratic f{R) functions are particularly useful 
to model inflation. 

The stability conditions are satisfied in the following 
manner respectively: 

• Conditions (1) and (3) give /' = /o|v^ > and 
/" = ./ojT^ > since R is positive. Condition (4) is 
not valid in this case. Condition (5) yields < ^ < 1 
at r = -2. 

• First two conditions hold /' = fovR""-^ > 0, /" = 
foi^ii^ - 1)R''^^ > whenever fo > 0, 1/ > I. Con- 
dition (3) is not valid here since A = and j^ 7^ 1. 
Condition (4) is invalid: it implies < (z^— 1) < 1 at 
r = -2 but ly > 1. 



The stability conditions f or f(R) gravity are (|] 
2O10l : lBamba et allboOStlApplebv et ail20ld ) 



vajf 



nclusion 



1. 



The stability conditions f or f(R) gravity are jiLuca fc Tsuiikawa 

In this paper we have considered the 4-dimensional 
f{R) theory of gravity. We have taken the spatially 
flat FRW Universe in the framework of f{R) grav- 
ity. Our aim was to investigate models which exhibit 
Noether symmetry with gauge term, where we have ig- 
nored the matter contribution in the action. On using 
the Lagrangian (7) in the definition of NGS (10) we 
have obtained a system of partial differential equations 
(14) - (20), which involve four unknown functions ^, 



/' > for R> Rq, where i?o > is the Ricci scalar 
at the present epoch. This is also required to avoid 
anti-gravity 

/" > for i? > i?o- It ensures consistency with lo- 
cal gravity tests, presence of the matter dominated 
epoch and the stability of cosmological perturba- 
tions. 



7], P and A, where each is a function of three vari- 
ables t, a and R. The integration of the equations 
(14) - (20) has yielded two different solutions ([26| and 
(|30p . Solving the system of equations, we obtained the 
gauge function to be zero. We then constructed the 
symmetry generators and the corresponding conserved 
quantities (also known as first integrals or constants 
of motion). We note that for the fractional power 
law, the algebra is three-dimensional and for the ar- 
bitrary power law, it is two-dimensional and thus the 
symmetry breaks by one. A previous study of ours 



(| Jamil et all 120111) discussing the /(i?)-tachyon model 
also resulted a zero gauge function. Thus we conjec- 
ture that the application of NGS to generic f{R) La- 
grangian results in zero gauge function. The symmetry 
generators give the time translational and scaling sym- 
metry of the theory. The stability analysis of the two 
forms of f{R) obtained here has been done in Section 
IV. We have shown that the two forms of f{R) given 
by ((26)) and ([30)) are stable and consistent with the lo- 
cal gravity tests by imposing some restrictions on the 
constants involved. Recently some black hole solutions 
have been studied in f( R) theory of gravity u si ng some 
known forms of .f(R) (Mvung et al gOllaU bl: iMvune 



20 lit iMazharimousavi fc Halilsovll2011 ). It would be 



of interest to study black hole solutions using the forms 
of f{R) obtained here. 
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